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Abstract 

The singularity structure of cosmological models whose matter 
content consists of a scalar field with arbitrary non-negative poten- 
tial is discussed. The special case of spatially flat FRW space-time 
is analysed in detail using a dynamical systems approach which may 
readily be generalised to more complicated space-times. It is shown 
that for a very large and natural class of models a simple and reg- 
ular past asymptotic structure exists. More specifically, there exists 
a family of solutions which is in continuous 1-1 correspondence with 
the exactly integrable massless scalar field cosmologies, this correspon- 
dence being realised by a unique asymptotic approximation. The set 
of solutions which do not fall into this class has measure zero. The 
significance of this result to the cosmological initial value problem is 
briefly discussed. 



1 Introduction 

The global structure of solutions of Einstein's equations, and in particular, 
the apparent existence of a cosmological space-time singularity has long been 
of interest to cosmologists. Considerable research effort during the 1960s and 
early 70s culminated in the detailed description of the cosmological singular- 
ity given by Belinskii, Khalatnikov and Lifschitz (BKL) ^, |^ and in the 
well known singularity theorems of Hawking, Penrose et al [^]. The latter 
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provide rigorous proof of the existence of singularities under very general cir- 
cumstances, and have long been regarded as demonstrating that singularities 
are inevitable in realistic cosmologies. 

In the absence of a complete theory of quantum gravity, the classical sin- 
gularity remains the only unambiguous means of describing the boundary of 
classical space-time and is therefore of importance in the formulation of a well 
posed initial value theory for cosmology. The cosmological initial conditions 
play a particularly prominent role in inflationary cosmology, which aims to 
provide a mechanism whereby the universe evolves from an arbitrary initial 
state toward the highly symmetric, isotropic universe we observe today [^]. 
At present, the question of whether or not inflation sucessfuUy achieves these 
goals remains unclear ^, largely, it is suggested, due to the difficulty in 
assigning meaningfull initial conditions in the absence of a global description 
of solution space. 

The global description of inflationary theories is complicated by the fact 
that a number of inflationary scenarios, such as chaotic inflation place the 
onset of inflation close to the boundary of the Plank era. On the other hand, 
inflation itself seems to arise from an essentialy classical mechanism and can 
be modelled by the evolution of a scalar field in a classical space-time back- 
ground. Both the classical and quantum mechanical limiting cases are there- 
fore of interest. It must be noted that the characteristic pre-infiationary time- 
scale for scalar field models will depend on the scalar field self-interaction 
term as well as some measure on the set of initial conditions, and there is no 
a priori reason why it should be of the order of the Plank time. 

The structure of the classical singularity requires careful re-examination 
in the infiationary context since the strong energy condition 0], which plays 
an important role in both the singularity theorems and the construction of 
the BKL singularities, can not be strictly applied in the early universe if 



infiation occurs 10 



In fact, it can easily be demonstrated by explicit example that there exist 
scalar field cosmologies which possess no singularity at all (see for example 
[pnp. Other well known solutions (ie, the power law infiation solutions ) 



possess singularities but no particle horizons |]T2| and it is possible, by pru- 
dent choice of the potential, to engineer exact Friedmann Robertson Walker 
(FRW) solutions with scale factor a{t) exhibiting almost any type of behavior 
one desires. 

In the present paper, we shall consider scalar field models generated by 
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an action functional of the following form: 



S[g,u. <!>]= J d^x^{R + + (1.1) 

where R is the Ricci scalar, g is the determinate of the metric g^i, and V is 
an arbitrary non-negative potential. (Note that we shall always work with 
units c = SnG = 1 and adopt the metric signature ( — h ++).) 

Although in what follows gfj_u is implicitely taken to be the physical metric 
on space-time, it should be kept in mind that effective scalar field theories of 
the form ( |1 . 1| ) arise naturally, via conformal transformations, in the context 



of generalised scalar tensor theories and higher order theories of gravity []T3 



ig, |T3|, |TB], 0, |TB|, ig, ^ 23). For such theories g^,^ will be equivalent to the 
physical metric up to a conformal factor determined by 0. 

Quantitative calculations of global behavior of the solutions generated 
by ( p..l| ) requires a detailed knowledge of the global form of the scalar field 
potential, a task that is fundamentally difficult since the energy density of the 
scalar field can be expected to diverge during the approach to a singularity. 
Interactions modeled by a simple power law, for example, would seem unlikely 
to remain valid over such a range. Furthermore, there are no known physical 
fields which have emerged as strong candidates for the "inflaton" field present 
in the early universe. It is by no means clear whether the scalar field driving 
infiation should be viewed as a fundamental physical field at all, or merely a 
convenient heuristic model which elucidates certain (qualitative) dynamical 
features of the early universe. The close relationship between scalar fields 
and conformal transformations, coupled with the likely significance of higher 
order gravitational lagrangians at high energies, is suggestive of the latter 
point of view. 

In accordance with the above remarks it shall be our approach to place 
strong emphasis on the structural features of scalar fields as a class of models, 
rather than the quantitative details of specific physical scenarios for inflation. 
We shall approach the problem from a dynamical point of view. If one insists 
that the energy density of the scalar field be non-negative everywhere, then 
it can be shown that the level hypersurfaces of the scalar field are space-like 
and (almost) always foliate space-time Therefore, a global cosmological 
time arises naturally and one can consider the evolution of the induced 3- 
metric and scalar field, together with appropriately defined time derivatives 
(or conjugate momenta), as an infinite dimensional dynamical system. 
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The easiest way to make further progress is by making use of the long 
wavelength approximation in which one neglects terms of second order in the 
spatial derivatives, thereby de-coupling the spatial and temporal dependence 
of the dynamical variables . The evolution at each spatial point is then 
governed by a finite set of ordinary differential equations equivalent to the 
evolution equations for spatially homogeneous space-times. The underlying 
physical assumption is that kinetic terms (ie, expansion and shear) diverge 
faster than spatial gradients and therefore dominate at very early times. 
For a more detailed discussion of long wavelength approximation schemes in 



inflation and cosmology in general see, for example, ^ ^ 

The long wavelength dynamics are characterised by the spatially homoge- 
neous scalar field cosmologies, the solutions of which reside in a 7-dimensional 
phase space. It turns out, however, that many of the interesting asymptotic 
features characteristic of scalar field cosmologies are already present in the 
relatively simple spatially flat Friedmann Robertson Walker space-time mod- 
els which, in addition to being an important subset of spatially homogeneous 
metrics in its own right, has the advantage that it can be represented by 
a two dimensional vector field. This allows one to visualise the flow and 
thereby gain considerable insight into the geometric structure of solution 
space, before tackling the higher dimensional problem. 

In particular, a number of authors have made fruitful investigations of 
inflationary structure in various FRW scalar field models using phase plane 
techniques (||2^, |l6l) and, importantly, this approach has proven to be 
readily generalisable to more complex space-times (eg, [0). Accordingly, the 
remainder of the paper shall be devoted to a detailed analysis of the general 
past asymptotic structure of spatially fiat FRW scalar field cosmologies from 
the geometric point of view. It is claimed that the techniques developed 
below can readily be generalised to more general space-times as is discussed 
in greater detail in P|. 

In Sections 2-4 we write down the field equations, assuming a general non- 
negative potential V, as a two dimensional dynamical system and investigate 
some of the elementary features of the flow. It turns out that the scalar field 
almost always diverges asymptotically in the past making it necessary to 
introduce some additional regularity conditions on V (which is unbounded 
in general) that allow us to compactify phase space and analyse the structure 
of the system at infinity, without the need for an explicit choice of potential. 
In Section 5 it is shown that under very general conditions there always exists 
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a generic class of solutions with an initial singularity which possesses a simple 
and elegant structure. These solutions are in continuous 1-1 correspondence 
with the exactly integrable massless scalar field cosmologies {V = 0), which 
provide a unique asymptotic approximation. In Section 6 these results are 
used to prove a global singularity theorem for FRW scalar field cosmologies 
and the cosmological significance is briefiy discussed in Section 7. 

Throughout the text major results are summarised as theorems. Refer- 
ence is made to a number of results from dynamical systems theory. Readers 
unfamiliar with dynamical systems techniques, and, in particular, the prop- 
erties of center manifolds are referred to the excellent introductory texts by 



Guckenheimer and Holmes p9[ and Wiggins which include comprehen- 
sive lists of major references. 

2 The Field Equations. 

Consider the FRW line element 

3 

ds^ = -dt^ + a{tYY.idx'f- (2-1) 

1=1 

We define the expansion K{t) as the function 

K = 3-, (2.2) 

a 

where the "dot" indicates differentiation with respect to time t. K may be 
interpreted physically as the rate of expansion of the spatial volume element. 

The evolution equations obtained from variation of (|1 . 1|) may be expressed 
as a dynamical system in the variables {K, 0, 0) as follows |^ : 



dt 

dK 
'dt 

subject to the constraint 



-K<j) - y (0) (2.3) 



= 3V(0) + ^02. (2.4) 
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The constraint equation actually makes ( |2.3| c) redundant since K can be 
expressed function of and 0. 

It can easily be verified that, for non- negative potentials, the set = is 
invariant under ( p.3|) with constraint ( p.4|) . Since no two orbits may intersect 
it follows that the sign of K is invarient. The physical expanding solutions 
of ( |2.3| ) therefore consist of trajectories in the 2-dimensional phase space 

n = {{K,(f),<i)) : K >0,K^ = 3V{(j)) + f^^} C R^. 
It shall be convenient to replace K and cj) with the new set of coordinates 

^ = ^ y = ^2^ (2-5) 

and introduce a new time coordinate 

T = \nv{t) + To (2.6) 

where f = is the spatial volume element and tq is some constant, r is 
well defined since v is strictly increasing. Furthermore, v goes to zero if and 
only if K goes to infinity 0] so ( |2.6|) implies that i^^ooasr^— oo. 
Differentiating ( p.6|) we find (recalling that K = v/v) 

dt dr ^ ^ 

In terms of these new coordinates the dynamical equations for the system 
become 



dx 
d^ 

dy_ 

dr 
dcf) 

dr 

subject to the constraint 



-y-Jlx'V'{ct>)+y' (2.8) 



+ 3xV(0) = 1. (2.9) 



If we imagine tracing a trajectory of the above system backwards in time 
then it may be seen from the above equations that both x and y remain 
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bounded. However, it is not possible to eliminate from the above equations 
since the constraint ( p.9|) is degenerate at x = and = 0. The analysis 
of the above system under the condition that remained bounded would 
be relatively straitforward, particularly since within some bounded region it 
would be appropriate to model V by a polynomial of sufficiently high order, 
greatly simplifying the problem. If diverges, on the other hand, then any 
such approximation schemes will almost certainly be inappropriate. It is 
therefore of considerable interest that almost always diverges as shall now 
be shown. 

Theorem 1 Assume that V is C^. Letp be a point in Q and let 0"{p) be the 



past orbit of p under (\2.d^ ). Then is almost always unbounded on O (p). 
Proof: 

Let p G n be such that is bounded on 0~{p). It is therefore clear that 
0^{p) is contained in a compact subset of (the closure of) Q. It is intu- 
itively obvious, and indeed follows from a basic result of dynamical systems 
theory (See [^], Proposition 1.1.14) that the trajectory must asymptotically 



approach some limit set a{p). It is easy to convince oneself that a{p) must 
itself be a union of orbits. Since x is monotonic, it must clearly be constant 
on a{p). ( |2.8| a) therefore implies that either y = or x = 0. The latter case 
would imply, by (|2.9|) , that = 1 on a{p) in which case it would follow 
from (p.8|c) that is unbounded, in contradiction to our previous assertion. 
Therefore, a{p) lies on the plane y = 0. 

Inspection of ( p.8| ) reveals that the only possible invariant sets lying on 
the line y = are equilibrium points (a;o,O,0o), where V'{xo) = and 
xl = |V^(0o) > 0. 

In order to complete the proof we need only to demonstrate that no 
such equilibrium point is a past asymptote to an open set of trajectories on 
Q. Note that the equilibrium points of (|2.8|) correspond to the equilibrium 
points of ( |2.^j| ), which may be expressed explicitely as a 2-dimensional (un- 
constrained) autonomous system by making the natural projection onto the 
plane 

X = 0, r = 0. (2.10) 



and eliminating K, using ( |2.4| ), to obtain 



X = Y 

/ o Xi (2.11) 

Y = -V'{X)- (3V{X) + IY^Yy. ^ ' 
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Let q = (Xo, Yq), be an equilibrium point of ( |2.11| ) satisfying 



Fo = V'{Xo) = V{Xo) = Vo>0 (2.12) 
We can classify the stability of equilibrium points by computing the eigen- 



values of the matrix of derivatives of the vector field ( p.3|) at q . We find 







The eigenvalues of Jg are: 



3Vo±j3Vo-AV"iXo: 



(2.13) 



(2.14) 



At least one of these eigenvalues must always be negative. It follows from 
the Center Manifold Theorem (See ||2^, Theorem 3.2.1) that there exists a 
local stable manifold intersecting q. That is, a 1 or 2-dimensional invarient 
submanifold of Q on which all solutions close to q exponentially approach q 
as t — > oo. The existance of a stable manifold of dimension n > implies 
that all solutions past asymtotic to q must lie on an unstable manifold or 
center manifold of dimension 2 — n < 2Q This completes the proof. □ 



It is of interest to dwell briefly on the physical meaning of the equilibrium 
points, above. Firstly observe that they correspond to maximally symmetric 
solutions of Einstein's equations. With Vq > 0, g represents a vacuum de 
Sitter solution with scale factor 




The sign of the eigenvalues depends on the sign of V"{Xq). If V"{Xq) is 
positive (Xo a local minimum of V) or negative (Xq a local maximum), then 

^If the second eigenvalue, above, is non-zero, then q is a hyperbohc fixed point and the 
result follows from the Hartman-Grobman linearisation Theorem ([^, Theorem 1.3.1 ) 
If there exists a zero eigenvalue, then the topology of the flow near q is non-trivial and 
is characterised by a unique local center manifold which acts as an exponential attractor 
in the forward time sense and therefore a repeller in the reverse time sense (See jjsij. 
Corollary 4.7). Briefly a center manifold is the non-linear analogue of the center subspace 
from the theory linear ODEs. It can be loosely defined as an invariant manifold to the 
fiow which intersects q and is tangent to the zero eigenvector at q. 
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q will be a stable sink or a saddle respectively. In the latter case exactly 
two solutions are past asymptotic to q. It is interesting to note that these 
solutions correspond physically to unstable de Sitter solutions, similar to 
those discussed by Barrow [|^. They asymptotically approach de Sitter 



space in the past and possess neither singularities or particle horizons. 

If V"{Xq) = 0, one of the eigenvalues is zero and the behavior near q will 
be characterized by a local center manifold {X, Y{X)), tangent to the X-axis 
at q, whose dynamics will be dependent on V. 



3 Conditions on the Potential 

Theorem |l| makes it clear that in order to investigate the generic past asymp- 
totic behavior of the system ( ^.31 ) we need look at the behavior near (p = ±oo. 
In order to do this we shall find it necessary to make some further restrictions 
on the behavior of as diverges. 

Firstly we shall restrict attention to potentials which are non-zero for 
large 0. In particular, this excludes the case where V has compact support. 
This will not lead to a significant loss of generality however, since equations 
( |2.3|) with V identically zero are exactly integrable and the solution is well 
known (we shall write it down explicitly later on). 

Non-zero potentials may be classified qualitatively according to their de- 
gree of "flatness" which we choose to characterise by the ratio V'/V, where 
V is the derivative of V. Intuitively we expect that the flatter potentials 
are more likely to exhibit inflation, for example. Indeed, the condition that 
V'/V be small plays an important roll in the slow roll approximation and is 
often stated as an explicit condition for inflation . 



Exponential type behavior would be characterised by V'/V tending to a 
non-zero constant (as — > ±oo) whilst V'/V tending to zero would indicate 
slower than exponential behavior, such as a polynomial. Extremely steep 
potentials, or those with oscillatory type behavior, would be without a well 
defined limit for V'/V as ^ oo. 

In the remainder of this section we focus on the behavior of V near = 
-I-CX3, an analogous discussion applies near = — oo: 

Definition 1 Let V : R —>■ R be a non-negative potential. Let there exist 
some 00 > for which V{(j)) > for all > 0o and some number N such 
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that the function Wy : [(po, oo) —>■ R, 



-N 



V{<P) 



satisfies 



lim Wv{(p) = 0. 



(3.1) 



Then we say that V is Well Behaved at Infinity (WBI) of exponential order 



It is important to point out that N may be 0, or even negative. Indeed 
the class of WBI functions of order is of particular interest, containing all 
non-negative polynomials and other 'fiattish' functions. 

A trivial consequence of the above definition is the following: 

Theorem 2 let V he WBI of exponential order N then, for all \ > N , 



Wc shall now define a procedure for classifying the smoothness of WBI 
functions at infinity. We simply say that a WBI function is of class k if there 
exists a sufficiently smooth coordinate transformation z = f{4>), which takes 
infinity to the origin, for which Wy is mapped to a C'' function in some 
neighbourhood of the origin. Furthermore, we should only admit coordinate 
transformations / which are useful from the point of view of dynamics, in 
that they induce a smooth map on the tangent bundle (ie smooth vector 
fields are mapped to smooth vector fields). Therefore the first derivative of 
/ with respect to 0, that is /', should also be expressible as a C*^ function of 
z. 

We take note of the following convention, which will be assumed hence- 
forth. If we have some coordinate transformation z = f{(f)) which maps a 
neighbourhood of infinity to a neighbourhood of the origin, then if is a 
function of 0, g is the function of z whose domain is the range of / plus the 
origin, which takes the values; 



N. 



lim e-^'^y(0) = 0. 



(3.2) 




z>0 
z^O 
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Definition 2 A C'^ potential V is class k WBI if it is WBI and if there 
exists 00 > and a coordinate transformation z = f{(j)) which maps the 
interval [0o, oo) onto (0, e], where e = /(0o) o,nd lim f = 0, with the following 

additional properties: 

i) f is C''^^ and strictly decreasing. 

a) the functions Wy {z) and f (z) are on the closed interval [0,e]. 

We designate the set of all class k WBI functions £^ . The third requirement 
in the above definition, that the first derivatives of Wy and /' vanish at the 
origin, is an additional condition on the transformation f which has been 
imposed for reasons of convenience that will become clear later. It will not 
effect the class of a potential since if the first derivatives of /' exist at the 
origin (for some coordinate transformation z) they can always be made zero 
by choosing the new coordinate z' = z^ . 

As illustrated by Table I, a diverse range of qualitative behavior is incor- 
porated within the framework of potentials. It can easily be verified that 
all non- negative polynomials and linear combinations of exponentials are . 
The very general inflationary potentials discussed by Barrow and Barrow 



and Parsons are also . Furthermore, sufficiently "smooth" distortions 



to the shape of the potential should not change the situation very much. 
For example, it can be shown that any potential for which Wy is will be 
provided that the second derivative of Wy vanishes faster than 0"^^+'^) 
as ^ oo, for some e > P]Q This is a sufficient but not a necessary 
condition, as can be readily verified for the quasi-exponential potential 

V{(j)) = e^-^'t' 

using the transformation z = l/(ln0)^. 

dp 

As a final note before proceeding, observe that the condition — — (0) = 

dz 

in Definitional may be expressed equivalently (using the chain rule and inverse 
function theorem ) as 

/" 

lim — = 0. 

</>-*oo /' 



f Conversely, if W^^- > 0{<j>~^) then Wy > 0(ln</>) -> oo as </> cx). 
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V{ct>) 




z = 


W^y (^) 


f (^) 




A 

n 


0" 




0^ 2 


2 


1 ^3 




P 












2e 










1 ^5 


In 




(01n0)-^ 


(ln0)-i 


1 

ze z 


2 

— ^ 




f 




20"! + (01n0)-i 


(In 0)^1 


(2 + z)e-^ 


2 

— ze z 



Table 1: Simple examples of WBI behavior at large 0. n and A are arbitrary 
constants. 

It follows that l/f is WBI of exponential order and hence by Theorem ^ 
that it is exponential dominated. This implies in turn that 1// is also expo- 
nential dominated. The function / must therefore obey the following condi- 
tion: for all e > 0. 

hm — = hm — = 0. (3.3) 

<p—>00 J 0— >00 J 



4 The Flow Near 6 = oo 



In this section we shall characterise the qualitative structure of the flow near 
= 00. As in the previous section, analogous results hold near = —oo. We 
shall assume that V G £\. Define the set C f2 to be the set of points in 
VL for which > e^^ where e is any positive constant, e is chosen sufficiently 
small so as to avoid any points where ^ = 0, thereby ensuring that Wy is 
well defined. We now make the coordinate transformation 

(x,?/,0) h-^ (x,?/, z), z = f{(t)) (4.1) 

on fie, where /(0) tends to as ^ oo and has been chosen so that condi- 
tions i)-iii) of Definition ^ are sastisfied, with k = 2. 

Substituting these new coordinates into (p.8|) and making use of the con- 



straint equation (2.9) to eliminate x we obtain the 2-D dynamical system 



-I = -iy + j=,iWv{z)+N)){l-y') 



dr 
dz 

dr 



f'iz)y 



(4.2) 
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We may identify with its projection into so that we have: 



a = {0<2</(e-i),-l<y<l} 

The variable x can now be treated as a function on defined by the con- 
straint equation which becomes. 

+ 3xW{z) = 1. (4.3) 

The directional derivative of x along the flow generated by (^.2|) may be 



obtained directly from Equation (2.8a). 

Since /' and Wy are at z = we may extend ( ^.2| ) onto the boundry 
of to obtain a system on the closure of fl^, that is Q^. 

Recall also, from Definition that /' and Wy vanish at the origin and 
are each second order or higher in z, and that /' is negative on Q^. 

As a simple example, consider the special case of the power law potential 

v{^) = -r- 

n 



where A and n are positive constants. Introducing the coordinate z = 
and substituting the expressions Wy and /' given in Table 1 into equa- 
tion ( ^.2| ) gives the system 

dr ^ ^ ^ ^ ^ ' 

1 3 

Tr = -t/'- 

Observe that it is perfectly regular everywhere on Vt^ (in fact, everywhere on 
i?2). 

Returning to the general schematic representation of Vt^ is shown 

in Fig. |1]. The boundary dVL^ is the union of the unphysical sets dVt^i, dfl^^ 
and dQ^4, as well as the physical set dQ^2 on which the observables K, 
and are finite, dfl^i, dQ^^ and dQ^4^ are each invariant manifolds of ( f4.2| ). 
dQ^2 is the only element of the boundary which is transverse to the flow. It 
is through this set that physical trajectories enter and leave Q^- 
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Figure 1: Schematic representation of the compact set fl^ showing fixed 
points (with N=0) and the sense of orbits on the boundary. 



By (|4.2| b) it can be seen that the sign of ^ is opposite to the sign of y. 



Thus trajectories may enter Cl^ only on the right half of the plane and leave 
only on the left half of the plane. 

The points (— ;^A^, 0), (1, 0) and (—1, 0) are the only equilibrium points, 
and shall be referred to as pn, P+ and p_ respectively. Note that the equilib- 
rium point Pn only occurs for potenials of exponential order |A^| < ^/Q. In 
what follows we shall always assume this to be the case. It turns out that 
this bifurcation represents a transition between infiationary models and non- 
infiating models as is well known from the structure of exponential potential 
models (eg, p8|). 

The structure of the fiow ( [1.2| ) close to the equilibrium points p± can 
most readily be deduced from the total derivative of the vector field (^, ^) 
evaluated at the respective equilibrium points: 




(4.5) 



The eigenvalues of J{p±) are 2(1 — -^) and 0, with corresponding eigenvectors 
iy = (1, 0) and = (0, 1) respectively. 
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It therefore follows from Fig. |^ that: dQ^i contains local unstable mani- 
folds for p+ and p- (ie, solutons close to p± having z = are exponentially 
repelled from the respective equilibrium point); dQ^^ is a center manifold of 

and; dQ^^ is a center manifold of Recall that by center manifold, we 
just mean an invariant manifold tangent to the zero-eigenvector. Because 
the linear part of ( |4.2| ) vanishes, the restricion of the flow to such sets is slow 
compared to the exponential decay/growth tangent to the y-axis. 

Consequently, all solutions in the neighbourhood of p+ or p_ will rapidly 
approach the respective center manifolds as they evolve backwards in time, 
and the stability of p+ and p- depends on the asymptotic behavior of the 
solutions on dQ^s and dQ^A respectively (0], Theorem 2.12). 

We may write down the solution of ( |4.2| ) on dQ^s and dQ^A by solving the 
differential equation 

| = ±TO\/|- (4-6) 
Recalling that z = f{(f)) this is just 




which yields 

0=^(±r + 0) (4.7) 
where = 0(0) is a positive constant. Thus 

z = /(v/|(0±r)). (4.8) 

The flow is directed towards the y-axis on 80,^4 and away from the y-axis on 
9f2e3 as shown in Fig. [l|. Thus, p+ is a saddle (Fig. ^ja) and p_ is a source. 
(Fig. |b). 

We now turn our attention to the remaining equilibrium point p^ = 
{—^, 0). Evaluating the derivative of the vector field at this point gives 

-"""'"("'^^ o)- (") 

The eigenvalues are — 1 + ^ and 0. This implies that there locally exists a 
1-dimensional center manifold Cat through p^, which is tangent to the 2;-axis 
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z 



z 




Figure 2: The stability of the equihbrium points and p_. 

at pn (and, in general, only at Pn)-, and a 1-dimensional stable manifold 
tangent to (and indeed a subset of) dVLei- In accordance with our previous 
comments on the properties of center manifolds it is clear that Cat is an 
exponential attractor on a sufficiently small neighbourhood of and it is 
intuitively obvious from the geometry, and indeed is rigorously true ([pi|]. 
Corollary 4.7), that any solutions past asymptotic io pn must lie on the 
center manifold; ie: 

Theorem 3 Assume V E S"^ with N'^ < 6. All solutions asymptotic in 
the past to p^ He on a unique orbit Cn (corresponding to a single physical 
cosmology) which is (locally) a center manifold of p^. 

For our purposes it is sufficient to observe that Theorem |^ implies that at 
most one cosmology can originate at pn and, therefore, that p- is the only 
generic source in dfl^- It turns out however, that the center manifold Cat is 
of fundamental physical interest in its own right, being intimately related to 
the inflationary dynamics of a particular model. corrosponds to the In- 
flationary Attractor emphasised by Liddle et al , and plays a role directly 



analogous to that played by de Sitter space in cosmologies with positive cos- 
mo logical constant. As such its significance goes well beyond the FRW "toy 
cosmologies" being discussed in this paper, since it acts as an attractor in 
the full space of Einstein's equations, and allows inflation to be characterised 
precisely in a very general context. 

Furthermore, the unique cosmology living on C^v accounts for virtually 



all known exact spatially flat FRW solutions (see p^. Table 1), including all 
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singularity free and horizon free solutions (excluding the trivial de Sitter cos- 
mologies). The most significant exception is, of course, the general solution 
for = (massless scalar field) which shall be discussed in further detail 
below. 

The fundamental dynamical role of Cat in scalar field cosmologies, and 
its relationship to inflation, shall be explored in more detail in a forthcoming 
paper. 

Before proceeding to the next section, we may summarise the results 
above by constructing a phase portrait: Fig. |] sketches the flow for the 
special case, = 0, Wv{z) > 0, illustrating the qualitative behavior and in 
particular the topological structure of the system near = oo. 




Figure 3: Schematic representation of the flow on 17^ for Wy > 0. 



5 The Space-Time Singularity 

In this section we investigate the behavior close to the generic source p_ in 
more detail, and demostrate that it corresponds to a space-time singularity. 
The past asymptotic behavior of solutions close to p_ will be approximated 
by the exact solution on the center manifold up to an exponential error 
term (this may be readily verified directly). Thus for r sufficiently large and 
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negative we may write 

y = -l + O(e^0 z = f (vl(0 - r)) + 0{e^n (5-1) 

For some positive constant 7. Note that evaluating the inverse f~^{z) in 
the above expression and using the mean value theorem and the expressions 
( |3.3D we find that 

0=y|(r + 0)+O(e^-) (5.2) 

for some positive /?, as we would expect from ( |4.7] ). 

Let us now express the asymptotic solution (|5.1|) wholly in terms of our 
observables. As ?/ — > — 1, equation ( |2.8| a) becomes 



^ = x(l + 0(e^O). (5.3) 
ar 

Integrating gives the first order expression 

X = XqC^ . (5.4) 

By dl) 



t = J xdr 

= xqc" + ti (5.5) 



which implies that t ^ tj as r — 00. Substituting ( ^.5|) into ( |5.4| ) and 
using ( p. 51) gives 

K={t-U)-\ (5.6) 

Note that i^T — >■ cxd as t — > tj, thus the source p_ corresponds to an initial 
space-time singularity. Substituting (|5l5| ) and ( ^.6|) into (|5.1|) (recalling the 
definition of y) we find 

fln^^^ = -yf(t-t,)-i (5.7) 

where c = xqc'^ is a positive constant. 

The asymptotic solution (|5.6|) , (|5.7|) is the exact solution of the system 



~M when V is identically zero; ie the massless scalar field. Thus there exists 
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a generic class of cosmologies which, in a sufficiently small neighbourhood of 
the singularity, behave, approximately, as though the matter content were a 
massless scalar field. 

In order to make the above statement more precise it is necessary to do 
two things. Firstly, we shall obtain an error estimate for equations (|5.6| ) and 
( |5.7| ) close to t = ti. Secondly, we shall demonstrate that close to t = ti there 
is a continuous 1-1 correspondence between massless scalar field cosmologies 
and cosmologies of the model V{(f)) which are past asymptotic to p-. We 
shall assume, without loss of generality, that U = 0. 

Theorem 4 Let V E have N"^ < 6. There exists a neighbourhood J\f{p-) 
of such that for all p G Af{p-.) the trajectory ipp approaches p- asymptot- 
ically in the past and the solution may be written: 

K = r^ + 0(ey(t)) 

cf) = -./|ln- + 0(tey(t)) (5.8) 

V a C 

^ = -\/lt'' + 0{ev{t)) 

where 

ey{t) = tV{-^\nt). (5.9) 
Note that from (|5.7| ), (|5.5| ) and Theorem ^ 



lim eV^""1/(;z(r)) = limt"1/(0(t)) = V a > ^N. (5.10) 



In particular this confirms that for A^^ < 6 the error terms above are domi- 
nated by the ffist order terms and, furthermore, for A^^ < |, the error terms 
tend to zero (uniformly, as we shall see below) as t — > 0. 
Proof of Theorem: 

We proceed by calculating a second order term for the function x. Any 
collection of higher order terms which are to be discarded will be denoted by 
a generic h. By higher order terms we mean, more precisely, that the ratio 
of h to all undiscarded terms tends to zero as r — *• —oo. We have 

1^,^. (5.11) 
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where we know that — 1 as r ^ — oo and x approaches the first order 
solution ( 
From 



y2 = l-3V{z)x'^ 

= l-3xlV{z)e'^^ + h. 



Thus, 



We will now show that 



dlnx 



dr 



1 - 3xlV{z)e^^ + h 



7 



+ h 



(5.12) 
(5.13) 

(5.14) 



where 7 = 2 — \ \N . Define the function T = Ve Since V has exponen- 



tial order N it follows from ( p.2| ) that 

V{z)e^^dT = J T{z)e"'^dT + h. 
Consider the indefinite integral 



/(r) = / T{z)e'^^dT 



Integrating by parts gives, 



/(r) = ^T{z)e^^\: 



7^3 



e^^T'{z)ydT 



TO 



h{r) + h{r) 



Consider the second term. Let 5{t) = sup^/<^ ^^1^^ then, recall- 

ing that y"^ < 1; 



V2 



e^^T{z)Wv {z)ydT 



7^3 

< 5{r)\I{r)\ 

< 6{T){\h{r)\ + \h{r)\) 
Sir) 



< 



l-5i 



T 



■\hir)\ 



(5.15) 
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for r sufficiently small. Letting tq go to — oo demonstrates ( p.l4| ), as required. 
We may now integrate (|5.13| ), using (|5.14|) and ( ^.10|) to obtain, 



X 



(5.16) 



Integrating again to obtain a second order expression for t we get 

t = xo{e^ - ^xlV{z)e^^) + h 

where the second term on the left hand side has been estimated in precisely 
the same way as J(r) above. This expression may be inverted, to second 
order, to give 



xoe 



t + 



1 



+ h. 



(5.17) 



substituting this into ( ^.161) and recalling that x = 1/K we find; 

K = t-^ + -V((f))t + h. 

7 



(5.18) 



In a similar manner we may obtain an error estimate for 0. Recall that 

d(f) rj 

which, may be integrated using (|5.12|) , to obtain 

t 



\n--^Vi<l>)t''^+h. 



(5.19) 



where ( ^.17| ) was used to replace r by t in the above expression. Similarly, 

' - iV(0)t) + h. (5.20) 



Clearly, the second term on the right hand side of ( |5.19| ) tends to as t — 0. 
This allows us to Taylor expand V about = • 



I In ^ to obtain. 



bWy 



I In. 



h 



where 6 is a constant. Substituting this equation into (|5.18|) , (|5.19|) and 
(|5.20|) completes the proof. □ 
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Theorem 5 Let V & £\ he of exponential order N"^ < 6. For sufficiently 
small neighbourhoods M{p-) of p-, the set of physical orbits intersecting 
N'{p-) is homeomorphic to the set of massless scalar field cosmologies; ie, 
For each c > there exists a unique cosmology satisfying (\5. d^) , and this 
corrospondence is continuous. 

Proof: 

Although the above theorem seems intuitively rather obvious, its proof is 
non-trivial due largely to the fact that precise details of the transformation 
between the trajectories of the system ([4 .21) and the physical solutions with 
respect to proper time t is not known. 

Firstly, recall that for each solution past asymptotic to p_ there exists 

c = xoe^ (5.21) 



where Xq and are as defined in ( |5.4| ) and ( |5.1| ). More precisely we may 
define 

Xo= lim e"^a;(r) 0= lim (r+A/|0(r)]. (5.22) 

Each of these constants depend on the point p at which we take r = but 
for a given solution ippir) a translation r i-^ r + tq (which moves p along its 
orbit ) leaves c invariant. This reflects the fact that each orbit corresponds 
to only one physical cosmology. 

In order to see this more clearly we make use of the fact that diverges 
and is strictly decreasing with r close to It follows that we may rep- 



resent any orbit as a graph {(f),y{(f))). Substituting ( |5.22|) into ( |5.21| ) and 
remembering that x is a function of state we have: 

c= lim e^'^x{(j),y{(f))). (5.23) 

Thus c depends only on the graph t/(0) (i.e. the orbit) and is independent of 
the initial point p. 

Let Mip-) be a small neighbourhood of p_ and let y{0] 0) be a continuous 
one-parameter family of orbits covering Af{p-). 6 is a bounded non- negative 
real parameter with the property that 62 > 61 implies y{d2](f)) > y{6i;(j)). 
Since two orbits may never cross, the above family of orbits will be well 
defined as long as A/'(p_) is sufficiently small that each orbit can be writen 
as a graph wrt (p. 6 may be viewed as an angular coordinate. 
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We shall now demonstrate that c is a continuous, 1-1 function of 6. This 
will immediately prove the first part of the theorem and the second part of 
the theorem (continuity) is then a trivial consequence of the continuity of 
the flow and the coordinate transformations ( |2.6| ), ( [4.1|) . 

Define the function 



C{e-A) = x{<P,y{e;<P))eV2'i> 

then c is, by definition, the limit of C. Taking the derivative of this expression 
with respect to and making use of (|2.8|a,b) we have 



(5.24) 

Also, from ([4.3| ) we may obtain an explicit expression for C 

c(^;0f = ^(i-y(^;0n- (5.25) 

It is immediately apparent from this expression that c is a continuous since 
the orbits y{6] (p) are uniformly continuous in 9. Also, since y is strictly 
increasing with 6 (for fixed 0), ( ^.251) implies that C{6; (p) is strictly increasing 
with 6. 

Let 6*2 > Oi and, for convenience, denote the graphs y and C along the 
respective orbits by yi, Ci and ?/2, C'2- Taking the derivative of the difference 
C2(0) - C'i(0) using (|5.24|) we obtain 

where we have used the fact that 1/2 ~ l/i > 0. Integrating from some 0o to 
00 we obtain 

In f Et^jKMM) > 0. 



from which it follows that 



c{e2) CM 

c(^i) Ci(0o) 
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since C{6; 0o) is strictly increasing. Thus c is strictly increasing with w and 
hence 1-1. 

This proves the theorem. □ 

Theorem |^ demonstrates that V is truly dynamically insignificant in the 
neighbourhood of the singularity p_ in that the family of solutions which 
asymptotically approach p_ are completely characterized by the solution 
space of the massless scalar field cosmological model. We shall discuss some 
cosmological consequences of this result in Section 7 but before doing so we 
shall make use of the above results to prove a global singularity theorem for 
scalar field cosmologies. 



6 A Global Singularity Theorem 

In order to complete the analysis of the global system it is, of course, neces- 
sary to also look at the behavior near (p = — oo. This can be done easily by 



observing that the equations (|2.3| ) are invariant under the transformation 

(0,0) ^-(0,0) V^U (6.1) 

where t/(0) = V{—(j)). Thus, for a particular potential V, the behavior of 
solutions of (|2.3| ) near = -co is equivalent (except of course for the sign of 
0) to the behavior of the system near = oo with the alternative potential 
U. Provided that U is in S^, the proceeding analysis on can be applied 
(for a suitable choice of e). In what follows, we shall denote the set of all 
potentials which are (class k) Well Behaved at both infinity and negative 
infinity as S^. We denote the exponential order of functions at oo and 
— oo by A^oo and A^_oo respectively. 

Note that if we define the set Q{xq) as the region of phase space for which 
X < Xq, then it follows from the fact that x is monotonic increasing that this 
set is equal to the union of its past orbits. The procedure outhned above for 
defining a coordinate system at — oo allows us to embed fl{xo) in a smooth 
compact three dimensional manifold Tj{xo) such that the vector field defined 
by (|2.8| ) may be smoothly extended onto S(xo). In order to do this we simply 
define a coordinate atlas as follows: Firstly, we define the interior of S to 
be the set {(x, 0) : < x < Xq, — 1 < |/ < 1} where the coordinates are as 
defined in (|2.5|). Clearly this set is bounded in x and y. On the open subset 
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of this set for which cf) > t for some sufficiently small e we define a second 



coordinate patch {x,y,z) according to ( [4.1|) . An analogous coordinate patch 
{x,y,z) is defined close to = — oo according to the procedure outlined 
in the proceeding paragraph. The construction of S(xo) is completed by 
attaching a boundary which is defined by taking the union of the planes 
X = 0,x = Xo,y = ±1, z = on the respective coordinate patches. 

Clearly S(xo) is compact. Indeed it has the topology of a big rectangular 
box so it could itself be embedded in (Although it is useful to keep this 
picture in mind, there is no need to carry it out explicitly.) The vector field 
( |2.8| ) can be smoothly extended onto the boundary of S(xo) such that E(xo) 
is the union of its past orbits. 

fl{xo) is a 2-dimensional hypersurface embedded in S(xo). The impor- 
tant feature to note is that fi{xo) approaches the (non-physical) boundary 
precisely along the intersection of the plane x = with the planes y = ±1 
and z = 0. We shall term this set, which has the topology of a rectangle, 
and is independent of xq, the unphysical boundary of Q, or dQ. Note, in 
particular, that it contains each of the equilibrium points p+,P- and pn at 

= ±00. 

Keeping all this in mind it is strait-forward to prove the following global 
singularity theorem. 

Theorem 6 Let V be a S"^ potential of exponential order < Q, then 

there exists a class of FRW cosmologies, containing almost all solutions of 
l \2.^ , which possess an initial space-time singularity. Furthermore, these cos- 
mologies are in continuous 1-1 correspondence with the exactly integrable set 
of massless scalar field cosmologies. This correspondance is realised asymp- 
totically according to the the following approximation 



±V|ln^ + 0(te^(t)) (6.2) 



where 



(Sketch) Proof: 



K = t-i + 0(e^(t)) 
= ±J^\n- + 0{ 

V cs g 

= ±v|r' + 0(e±(t)) 



e±(t) = tV{±J^\nt) (6.3) 
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Clearly, it is sufficient to show that almost all solutions are asymptotic 
in the past to the equilibrium point p- (at either oo or —oo). Also, since x 
is monotonic, it is sufficient to consider solutions on fl{xo) C S(xo) where 
Xq is arbitrary. Since E(xo) is compact and contains its past orbits, each of 
its points p must contain an a-limit set, a{p), which is itself a closed union 
of orbits. In particular, for points on the physical space Q{xq), Theorem |l| 
implies that ap must almost always contain at least one point on the sets 
z = {(f) = ±00). But, from the discussion proceeding the theorem, any 
point z = which is a limit-point of a physical trajectory must be part 
of the unphysical boundary dQ and must therefore have x = 0. Since x 
is monotonic increasing, the entire of a{p) must be contained in the plane 
X = 0, or more precisely, on the 1-dimensional "rectangle" dQ. It can easily 
be shown that the only conceivable generic limit sets on this rectangle are the 
two equilibrium points corresponding to p_ at ±00 respectively (By Theorem 
^ we know that Pat is not a generic source, neither can the rectangle be a 
cycle since the sense of orbits reverses at Pn)- 

This completes the proof. □ 



7 Discussion 

Perhaps the most interesting aspect of Theorem ^ is the indication it pro- 
vides that scalar field cosmologies have a simple and quite regular asymptotic 
structure, and that this structure is independent of the exact details of the 
potential. This allows one to impose initial conditions and physically moti- 
vated boundary conditions in an unambiguous manner without the need to 
make an arbitrary choice of initial spacelike hypersurface. Rather, one takes 
the boundary as being the singularity itself and assigns boundary conditions 
according to the unique asymptotic state of the solution. This corresponds 
to choosing a particular massless scalar field solution which, in the case of 
the very simple FRW models above, is equivalent to choosing a parameter 
value c (together with a choice of sign). 

For example, if one makes the assumption that all values of c are a priori 
equally likely, then it is possible to unambiguously assess the probability of 
inflation in a particular model by simply comparing the size of the line seg- 
ments in c-space which inflate and don't inflate, respectively. Alternatively, 
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one could carry out quantum gravitational calculations on the massless scalar 
field cosmologies to obtain a probability distribution P{c) on the set of ini- 
tial conditions. In view of the technical difficulties associated with quantising 
scalar field cosmologies with non-zero potentials, whose classical solutions are 
not generally integrable, the appeal of the above procedure is self evident. 

FRW space-time does not provide a particularly realistic model of the 
gravitational field close to a singularity due to the tendency of gravitational 
perturbations, in particular shear anisotropy, to diverge as the singularity is 
approached. The significance of Theorem ^ rests with the relative ease with 
which it may be generalised to more general space-time behavior (and more 
complicated matter fields). Investigations carried out so far have provided 
a fairly strong indication that the dynamical insignificance of V close to 
the singularity is indeed a characteristic feature of scalar field cosmologies 
in general. A direct analogue of Theorem ^ has been proven for Bianchi 
type I space-time and a weaker result has been demonstrated for Class A 
spatially homogeneous space-times, which include Bianchi type VII and IX 
- the highest dimensional spatially homogeneous cosmologies; ie, it can be 
shown that almost all solutions can be approximated by \^ = solutions 
asymptotically in the past, under the condition that a bound is placed on 
positive spatial curvature at some time to HI- 

In fact, the presence of positive spatial curvature introduces significant 
additional structure into solution space, due to the fact that the time deriva- 
tive of the expansion K need no longer be monotonic (as is well known from 
standard cosmology). For scalar fields this leads to the existence of singu- 
larity free "bouncing" solutions for which the spatial volume element has a 
non-zero minimum, even in relatively simple FRW models [^]. For the case 
of the massive scalar field {V = l/2m?(fp') these solutions have been studied 
by Page |^ and, more recently, by a number of authors |38| who have 
shown that they comprise a rather complex fractal subset of solution space. 
Although intriuging, this structure has measure zero in solution space and is 
therefore non-generic. In models for which V has a non-zero minimum, on 
the other hand, it appears that bouncing solutions are generic, not only in 
FRW models but in the more general Bianchi Type IX solution space. The 
possible physical significance of such solutions is yet to be looked at in de- 
tail. Nevertheless, as is the case with no-hair theorems for inflation p9| , it 
seems possible to eliminate cosmological bounces by placing a bound on neg- 
ative curvature. This is suggestive of a close relationship between classes of 
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past asymptotic structure and the occurrence of inflation. The past asymp- 
totic structure of the non-bouncing solutions appears to conform with that 
of negative and zero curvature models, in that neither the curvature, nor the 
scalar field potential are dynamically significant asymptotically in the past. 

Given that the long wavelength solutions of scalar field cosmologies asymp- 
tote towards the solutions of the massless scalar field {V = 0) it is significant 
that such a field is equivalent to a perface fluid with equation of state p = p; 
the so called stiff fluid. This extreme case among classical linear-equation- 
of-state perfect fluids differs from other such models in that its initial sin- 
gularity is believed to have a velocity dominated semi-Kasner like structure 
rather than the more characteristic, and much more complex, oscillatory 



Belinskii Lifshitz Khalatnikov singularity |0, ^|. It is worth noting that the 
semi-Kasner like singularity structure admits the long wavelength approxi- 
mation close to the singularity wheras the BKL singularity does not. In both 
cases however, the procedure for constructing the generic singularity is es- 
sentially to take the general spatially homogeneous solutions (with respect to 
an orthonormal frame, say) and make the parameters dependent on spatial 
coordinates. Thus, as far as we know, the asymptotic behavior of general 
cosmological perfect fluid solutions is very closely related to, and is indeed 
characterised by, the solution space of spatially homogeneous space-times. 

More work is required before we can state with confidence that this is 
also true of scalar field cosmologies but it seems probable that it will at least 
be true of a generic family of solutions, provided one is careful to insist that 
the scalar field has non-negative energy (ie; that the weak energy condition 
is satisfied). This condition may be justified on heuristic grounds, and has 
the effect of imposing significant constraints on the degree of inhomogene- 
ity of the scalar field. The author would further speculate, that the weak 
energy condition is also likely to be a necessary condition in ensuring that 
significant proportion of solutions undergo inflation in accordance with the 
already stated view that there exists a close relationship between singularity 
structure and inflation. 
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